Resonances in one and two rows of triangular Josephson junction 

cells 



We present an experimental and analytical study of resonances in the 
current-voltage characteristics of single and double row triangular Josephson 
arrays. The magnetic field dependences of the voltage positions of the reso- 
nances have been measured for various cell inductances and junction critical 
currents. In double row arrays, we have observed a peculiar resonance, whose 
voltage position decreases with magnetic field. We derive the spectrum of 
linear electromagnetic waves propagating in the arrays. In the double row 
array, the spectrum consists of many branches that differ by the behavior of 
the Josephson junctions transverse to the bias current direction. The mea- 
sured magnetic field dependence of the resonance voltages is mapped to the 
linear mode spectrum and good agreement between experiments and model 
is found. 
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I. INTRODUCTION 



The driven Josephson junction ladders and arrays have recently attracted a wide interest 
due to many fascinating phenomena observed in such systems. They show vortex propaga- 
tion various current- voltage resonances nonlinear dynamic localized modes 
and may also be practically useful as high-frequency oscillators 

A peculiarity of quasi-two-dimensional arrays, i.e. Josephson ladders, and two- 
dimensional arrays is the presence of Josephson junctions in both the longitudinal and 
transverse directions to the bias current [see sketches in Figs.|l|(a)-(c)]. In these systems, 
the Josephson phase dynamics is more complex than in the rather well studied case of one- 
dimensional parallel arrays due to the higher degree of freedom appearing from the 
presence of the junctions transverse to the bias. For example, the spectrum uj{q) of the linear 
electromagnetic waves (EWs) propagating in the ladders [Fig.|I|(a)] contains two branches 
1^,^]^. Similar behavior has been extensively studied in vertical stacks of long Josephson 
junctions [see Ref. [^ and references therein], where the inductive coupling between the 
junctions leads to the splitting of the dispersion relation. The splitting has been found also 
in coupled one-dimensional parallel arrays |T^. As we turn to the case of two-dimensional ar- 
rays, the number of branches in the spectrum uj{q) increases. This feature has been recently 
confirmed by a systematic numerical study on dynamical states in underdamped arrays [ITT 



It is well known that, in the presence of magnetic field, the nonlinear interaction between 
the Josephson current wave and the excited EWs leads to resonant steps in the current- 
voltage characteristics of extended Josephson systems |]I|,|13|,|l5 . The measured magnetic 
field dependence of the voltages of the resonant steps can be mapped to the spectrum u{q) 
by the relationship: 

K„ = M«). , = / . (1) 

where / = $ext/$o is the magnetic flux threading the cell normalized to the magnetic 
flux quantum, and it is referred to as frustration parameter. The spectrum uj{q) for the 
Josephson ladder with four junctions per elementary cell [Fig.|l](a)] has been calculated in 
Ref. 101 in the linear approximation, and a good agreement between the measured magnetic 
field dependence of the step voltages and Eqs. (|l]) was found. This system has a remarkable 
property that at / = 0.5 an "out- of -phase" state of the currents is formed. As the 
amplitude of the resonant step reaches a maximum, the direction of the mesh currents 
flowing in adjacent cells alternate from one cell to another, and the so-called " checkerboard " 
structure of the alternating (ac) currents appears @J^. The ac currents compensate each 
other on the horizontal junctions and the total alternating voltage along the ladder is close 
to zero. This property is an obvious drawback if one wants to use the ladder as an oscillator. 
To overcome it, Yukon and Lin proposed |0| to use a triangular arrangement of Josephson 
junctions along the ladder [Figs. |l](b)-(c)]. In this case, at / = 0.5 an " in-phase" state 
with all the ac currents flowing in the same direction can be realized for the junctions 
transverse to the dc bias current. So far, the detailed study of the EWs dispersion in 
triangular Josephson ladders and arrays has not been carried out, as well as the study of 
the dependence of the resonant steps on the magnetic field. Thus, in this paper we present 
measurements of the resonant steps in the current-voltage (/-V") characteristics of single 
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and double row triangular Josephson arrays. For both systems, we calculate the spectrum 
of linear EWs and map it on the magnetic field dependence of the step voltages. The paper 
is arranged as follows: in Sec. II our experimental observations are presented. In Sec. Ill, 
we introduce the model based on a set of equations, and derive the spectrum of linear EWs. 
Finally, Sec. IV contains a discussion of the obtained results and conclusions. 



II. EXPERIMENTS WITH TRIANGULAR JOSEPHSON JUNCTION LADDERS 

AND ARRAYS 

In order to investigate the dynamical behavior of triangular arrays in the presence of an 
external magnetic field, we have chosen linear arrays made of one row and two rows, as shown 
in Figs, ^(b)-(c). Arrays of various critical current densities and sizes of the elementary cells 
have been measured. The arrays were made of Nb/Al-AlO^^/Nb underdamped Josephson 
tunnel junctions, arranged in a triangular lattice. The junction area is designed to be 
9/im^. The studied critical current densities are jc ~ 50A/cm^ and jc ~ 1000 A/cm^. 
The junction capacitance C is about 300 pF and 450 pF, respectively. The bias current If, 
is uniformly injected in each node of the array via on-chip resistors and extracted as shown 
by arrows in the sketches. The voltage is measured in the direction along the bias, across 
each individual row. The measurements have been performed using the acquisition software 
Ref. [|l5l. Experiments have been performed in the presence of a magnetic field applied 
perpendicular to the cell plane. Following the standard notation, we express the magnetic 
field in terms of the frustration /. The self-inductance of one cell can be roughly estimated 
as 



L = 1.25/ioVA , (2) 

where /iq is the free space permeability, and A is the cell hole size. The value of L is needed 
for calculating the parameter (3^ = 27rLJc/$0) where Ic is the junction critical current. The 
value given by Eq. ^ underestimates the actual inductance of the cell, as it is an asymptotic 
value valid only when the width of the superconducting electrodes forming the cell is larger 



than the hole size [Ig]. On the contrary, our arrays are in the opposite limit. We have 
compared the calculated value (Eq. 0) with the value L obtained from the magnetic field 
dependence of the critical current of a single cell with two-junctions, i.e. a SQUID. The 
ratio of the minimum to the maximum critical currents of the SQUID refiects the value 
of its parameter |]14| and, therefore, the value of the cell inductance [0. The relation 



between the two values was found to be L ^ 0.7 L. However, in order to make systematic 
comparison of all studied arrays, the Pl values reported here refer to the cell inductance 
estimated from Eq. ^ To get various values of (3^, we have used samples either with different 
critical current densities and same cell size, or have varied the cell size A from 126 /im^ to 
240 /im^ for similar arrays located on one chip. Thus, the parameter varied between 
0.5 to 5, at the temperature T =4.2K. The investigated arrays and their parameters are 
summarized in Table I. 

The measured I-V curves of the triangular ladders show well defined resonances, grouped 
in two different voltage regions. Fig. ^(a) shows an enlargement of the I-V curve of a ten 
cell row with ~ 0.5 at / = 0.3. Here Si and S2 denote the upper and lower voltage 
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resonances, respectively. The steps of type 5*1 and 5*2 are stable in approximately the same 
range of frustration. Changes of / induce a periodic modulation of the voltages Vi^2 of the 
steps. In particular, Si and 5*2 approach the maximum voltage at / = 0.5 (and other half- 
integer values), and tend to the minimum value of voltage at integer /. FigureH(b) shows the 
^1,2 (/) dependences in the range < / < 1. The step S2 has a typical resonant behavior. 
It clearly shows distinctly different voltages (underlined by the dashed lines of Fig.|^(b)) 
that correspond to the minima of the differential resistance on the step. However, we have 
noticed that the resonant regime of the step S2 is realized only in the case of low (3l, i.e. for 
small discreteness. In arrays with either larger 1^ or larger L (both these quantities lead to 
an increase of Pl) the resonant behavior disappears and the step S2 changes continuously 
with /. 

The voltage and the differential resistance of the step Si have almost a continuous de- 
pendence on / for all the investigated (3^ parameters. In general, the differential resistance 
of the step Si is higher than that of S2. This feature is enhanced in the samples with large 
I3l, as Fig.|^ shows for the case of a single row array of 12 cells and (3^ ~ 5. At a constant 
bias current chosen on the step Si, by changing / we could continuously tune the voltage 
across the array, along the periodic pattern shown in the figure inset. 

The maximum voltage of the steps 5*1 and S2, i.e. the voltage value approached at 
/ = 0.5, depends on the cell inductance. We have compared the maximum voltages of three 
triangular ladders (samples #3, #4, #5) designed with different cell sizes {A = 126 /im^, 
180 urn^ and 240 /im^, respectively) and same junction area, and fabricated with the same 
critical current density {jc ~ 50A/cm^). The data reported in Fig. | show a tendency 
for the maximum voltages of both Si and 5*2 to decrease as the (3l (cell size) increases. 
This effect seems to be more pronounced for the step 5*2 than for the step Si. Increasing 
temperature also leads to a reduction of the maximum step voltages, similarly to the results 
reported in Ref. 

In experiments with two row arrays [Fig.|I|(c)], we have used a common bias current 
for both rows, and have independently measured the voltage across each of them. Figure]^ 
shows the typical I-V characteristics of an array made of 12(cells) x 2(rows), sample #6. 
Both rows (denoted as A and B) have very similar I-V curves. At the same critical current, 
the rows switch simultaneously to the finite voltage state, with the same voltage. Similar to 
the case of single row arrays, in the presence of frustration we distinguish the steps Si and 
S2. In Fig.|^ these steps are plotted at / = 0.4. No hysteresis is observed in this case. The 
steps are found to be periodically modulated by the magnetic field, in a similar way as it 
is described above for the single row array. Moreover, their maximum voltages approached 
at / = 0.5 are similar to that observed in the single row array with the same discreteness 
parameter (sample #2), i.e. with the same cell size {A = 160 fim^) and critical current 
density {jc = 1050 A/cm^). The new feature observed in the two row arrays is the presence 
of a third step, marked as S3 in Fig.^ The peculiarity of this resonance is that it is stable 
at frustration values corresponding to an integer number of fluxons per cell (/ = 0, ±1, ...). 
At these values of /, an increase of lb above the array critical current induces both rows 
to jump simultaneously to the state S3. Further increasing of causes the transition from 
S3 to the McCumber branch. The I-V curve is hysteretic and a decrease of h along S3 
eventually reaches an instability point with a certain retrapping current, at which the array 
returns to the zero voltage state. Sometimes, by decreasing /(,, at the instability point the 
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two rows split and while one row goes to = 0, the other row undergoes a transition to 
higher voltages and then to = 0, with a slightly lower retrapping current (this behavior 
is shown in Fig.^. 

The dependences of the maximum voltages of the steps Si, S2 and Sg on / are shown 
in Fig.^. In contrast to Si and S2, the step S3 moves to lower voltages when approaching 
/ = 0.5. We note that the voltage modulation is more pronounced for the steps Si and S2 
than for S3. Moreover, increasing the temperature above 4.2 K causes only a slight reduction 
of the step voltage. Above / ~ 0.2 the step S3 disappears and the steps Si and S2 become 
stable. 



III. MODEL AND SPECTRUM OF LINEAR ELECTROMAGNETIC WAVES 

The derivation of the equations of motion for a triangular array is done in a way similar 
to the case of ID-arrays (2 junctions per cell) |Q and ladders (four junctions per cell) 0,|T2 



The junctions are described by the RCSJ model |]T4[. We neglect mutual inductances and 
consider only the self inductances of the cells. We denote with (pij and (Pi+ij, respectively, 
the superconducting phase differences (Josephson phases) across the vertical junctions of 
the cell i in the row j {j = 1, 2) and with ipi^k the phase differences across the horizontal 
junction of the cell i in the line k {k = 1,2,3). For the one-row case the indexes j and 
k are unnecessary and, therefore, omitted. With $j and $ext we denote the induced and 
the applied magnetic flux, respectively. First, we derive the equations of motion for the 
"triangular" Josephson ladder of Fig. 1(b). We recall the fiuxoid quantization in the cell i: 

(fi+i -(fi + iJi = — — . (3) 

Due to nonzero cell inductance L, $j and $cxt are related by 

= <l>ext + Lli , (4) 
where is the mesh current in the cell i. By making use of Kirchhoff 's current law, we get 

CV^' + ^ + hsmif, = h-U + h.i , (5) 

where Vj" is the voltage across the vertical junction i. All junctions have capacitance C, 
resistance R and critical current Ic- Since the current flowing through the horizontal branch 
is the mesh current /j, for the horizontal junctions we obtain 

CV^ + ^ + IcSiniji = h . (6) 
Finally, the equations of motion for the vertical and horizontal junctions of the row read 

as: 

LPi + a^pi + sin Lfi = 

PL 

1 ^Tr^ 
i)i + aipi + sinipi = -riifi - ipi+i - ipi) - —— . (7) 

PL PL 



5 



Here, the time unit is tUp ^ = \JhC l{2el^, the inverse plasma frequency. The parameter 
a = Xj \f^c determines the damping fl^ of the junctions and the parameter /^^ defines the 
discreteness of the array; / is the frustration defined above, and 7 = /f,//c is the normahzed 
bias current. 

To derive the spectrum of EWs, we assume a whirhng solution along the vertical junctions 
and oscillations with a small amplitude for the horizontal junctions. Moreover, the phase 
of the vertical junctions increases from cell to cell due to the presence of frustration. Thus, 
the solutions of Eqs.|^ can be written as: 

Lp^ = ujt^ 2-Kfn + LP e^('^*+2^'?") 

^„ = ^ e^('^*+2"^") , (8) 

where uj and q are, respectively, the angular frequency and the wave number of the EW in the 
array. In the limit of small amplitudes (/? and ^Z', we obtain the spectrum of electromagnetic 
wave propagating along the array. This spectrum consists of two branches a;+(g) and a;_(g) 
given by: 



cu± = uj^^T±4W^, (9) 

where = 1/2 + 1/(2/3^) + (2//5l) sin2(7rg), and Q = {A/ /3l) sin'^{7cq). The two modes 
are plotted in Fig.^. This dispersion relation differs from that derived in Ref. for the 
case of "square" Josephson ladders [four instead than three small junctions per elementary 
cell. Fig. 0(a)] only by a constant factor. The horizontal junctions play an essential role 
in the array dynamics, and lead to two linear resonances in the dispersion relation. If, 
instead of the horizontal junction, simply a superconducting link is placed, that is the case 
of well-known ID parallel array, only one linear mode exists . 

Thus, what should we expect when two rows of cells join together in a 2D array? In order 
to describe the spectrum of linear modes in the two-row array, we use the time dependent 
Josephson phases of vertical v^ij(t) and ipi+ij{t) and horizontal ipi,k{t) junctions. Similar 
to the one row case, we derive the set of equations by means of the Kirchhoff's current law 
and the fluxoid quantization. Assuming a whirling solution along the vertical junctions and 
oscillations with a small amplitude for the horizontal junctions, we obtain the spectrum of 
linear modes from the system of seven linear equations: 

- /?Lt^Vi,i = V5i+i,i(l + e"*'') - 2v9j,i - - V'_i,ie"*'' 

-/5Lu;Vi,2 = V'i+i,2(l + e"*^) - 2ipi^2 + V^o,i+i + V'l.^e"*'' 
-/5Lu;Vi+i,2 = V5i,2(l + e*^) - 2v2i+i,2 - V^o,i+i - 



(10) 



By solving the system of Eqs. (p!OD, we obtain seven branches in the oj{q) dependence. Three 
branches are determined by the equation: 
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W| = . (11) 



4 ^ 29 



These branches correspond to a ribbon state 0, as the Josephson junctions of the middle 
row are not active [ipo = 0). The lower ribbon branches displays an uj{q) dependence which 
increases as the wave vector increases. As observed in the experiments, the upper branch 
displays a different behaviour with respect to the lower branches, i.e. uj increases as the 
wave vector decreases (alternatively, the voltage increases as the frustration decreases). 

The other four solutions of the system of Eqs. ( p!0D involve oscillations of the junctions 
in the middle row and, therefore, correspond to a checkerboard state (at / = 0.5). These 
four branches are determined by the equation: 

[(-..^ + |-)/3l(l-c.2 + -l)-2]x 

PL PL 

X [-2a;2 + (-^ - u^) (1 - uj^)] = 2(1 - uj'^f cos^ | . (12) 

PL 2 

As in the ribbon case, the upper branch of the checkerboard state has a different depen- 
dence on / with respect to the lower checkerboard branches. All the linear modes calculated 
for the double row array are plotted in Fig. As one can see, the qualitative behaviour of 
these linear modes is in a good agreement with the experimental data reported in Fig.|^. 
Thus, the linear approximation used for modeling the dynamic behavior of the frustrated 
arrays is sufficient to quantitatively explain the most essential features observed in the ex- 
periments. 



IV. DISCUSSION AND CONCLUSIONS 

In the previous section we derived the spectrum uj{q) of EWs for triangular Josephson 
ladders and for double row arrays. In the Josephson ladder case [Fig.|l|(b)], we obtain two 
branches uj±{q) [Fig.|^. The excitation of the EWs that account for these branches leads to 
resonant steps in the I-V curves [Fig. 0. Although the current amplitude of the resonant 
steps depends on the array parameters in a rather complicated way, the voltage positions 
of the steps can be mapped to the spectrum of linear EWs by making use of the Eq.(0). 
Indeed, we obtain a good agreement between the magnetic field (frustration /) dependence 
of the step voltages [Fig.|^(b)] and the calculated spectrum uj± (Fig. ^. The dependence 
of the maximum voltages on the discreteness parameter (5^ shown in Fig.^, displays also a 
good agreement with the theoretical prediction. Thus, the limiting operation voltage (i.e. 
frequency) for each mode can be controlled by the geometrical inductance and/or by the 
critical current density. 

The calculated spectrum of EWs for the double row triangular Josephson array contains 
seven branches (Fig.^. Three of them correspond to the ribbon state, while the others 
are due to the checkerboard state. However, the experimental data show that only three 
branches could be excited in the Josephson arrays and, correspondingly, only three reso- 
nances appear in the current-voltage characteristics. The reason of this might be ascribed 
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to internal instabilities of the states. The observed resonances 5*1, 5*2 and S3 (Fig. H) can 
be mapped to the spectrum of EWs (Fig.|^). Moreover, the observed peculiar resonance 
(5*3) appears at the small values of frustration /, and its voltage decreases with /. Two 
branches with the similar behaviour are also found in the spectrum uj{q) of EWs (Fig. 
one corresponding to the checkerboard and the other to the ribbon state. Unfortunately, 
measurements of I-V curves do not allow to distinguish between these branches and, corre- 
spondingly, between the ribbon and checkerboard states. A distinction can be, in principle, 
done by detecting radiation from the array due to properly coupled ac voltage. 

In conclusion, we have studied the dynamical states of triangular arrays of Josephson 
junctions in the presence of a magnetic field. As expected, the number of exhibited reso- 
nances increases with the number of degrees of freedom of the system. We found that in 
one-row arrays there are two states, and that in two-row arrays three states are observed. 
The proposed analytical model allows to obtain the spectrum of linear modes, both for sin- 
gle and double row triangular arrays. The voltage position of the observed resonances is 
mapped to this spectrum and good agreement between experiments and theory is found. 
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TABLES 

TABLE L The Arrays Parameters at T =4.2K. 



Sample" 


#1 


#2 


#3 


#4 


#5 


#6 


No. of rows 


1 


1 


1 


1 


1 


2 


No. of cells per row 


10 


12 


10 


10 


10 


12 


Cell area (/^m^) 


126 


160 


126 


180 


240 


160 


jc (A/cm2) 


30 


1050 


50 


50 


50 


1050 


Pl 


0.5 


5 


1.6 


1.9 


2.2 


5 


"All the studied samples 


have 3 junctions per cell. 
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FIGURES 



FIG. 1. Sketches of Josephson ladders with (a) four and (b) three Josephson junctions per ceh, 
and of a two-row array (c); ipij and denote the Josephson phases of the vertical junctions of 

cell i in row j; ipi^k denotes the Josephson phases of the horizontal junction of cell i in line k. The 
magnetic field is applied perpendicular to the plane of the cells and induces a mesh current 7^. 

FIG. 2. (a) I-V characteristic of a one-row array at / = 0.3; (b) step voltage dependences 
vs. / of Si and S2 (sample #1); dashed lines represent the well-defined resonances of the step 82- 
Here Pl =0.5 and T = 4.2 K. 



FIG. 3. Sequence of Si steps recorded while slowly changing / in the range 0-0.5. The single 
row array has 12 cells and a ~ 5 (sample #2). Straight arrows indicate the transition from the 
step to the McCumber branch. Inset: the voltage across the array at the bias point ~ 32 /xA is 
smoothly tuned from to 65 /xV as / changes. 



FIG. 4. Filled symbols: Maximum voltage of the steps Si and S2 measured at / = 0.5 in 
arrays having three different values of cell inductance L (samples #3, #4, #5). These arrays have 
10 cells, jc ~ 50A/cm^ and Up 7 GHz. Solid line is the theoretical prediction (Eq. 9) for the 
upper resonance, for q = 0.5. 



FIG. 5. I-V characteristics for the two row array showing the steps Si and S2 at f = 0.4 
(open and filled circles), and the extra step 5*5 at / = 0.1 (open and filled squares). Data refer 
to the voltages measured independently across the individual rows, noted as A and B. Straight 
arrows indicate the hysteretic path. Sample #6, (3l = 5. 



FIG. 6. Measured dependences of the maximum voltages of the steps Si , S'g and S3 as a 
function of frustration for the two row array (sample #6). Data refer to the voltage measured 
across one row (A), another row (B) exhibited simultaneously the same voltage state. 

FIG. 7. Linear modes a;+ and u- in one-row array with three junctions per cell (continuous 
line) and four junctions per cell (dashed line). /3l = 1. 

FIG. 8. Calculated spectrum for the linear modes in the two-row array and sketches of the ac 
mesh current distribution at / = 0.5. 
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